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Abstract 

Using a high-order accuracy finite-difference time-domain algorithm, the acoustic scattering 
from a flat-plate cascade is computed. Keeping the grid and time step fixed, the effect of four dif- 
ferent boundary conditions on the accuracy and stability of the computed solution is compared. 

Introduction 

Despite recent efforts to introduce quieter aircraft, it is expected that noise will increase as a 
barrier to air transportation 1 ' 3 . For high-by-pass jet engines, fan noise becomes comparable to jet 
noise 4 . With current emphasize placed on jet noise reduction, fan noise could become the domi- 
nant noise source in the near future. Fan noise can be classified into three main sources 5 : inlet 

boundary layer or inflow distortions interacting with the fan; self noise from the fan; and fan 
wakes interacting with stators or struts. Both tone noise and broadband noise can be generated. 

The gust-cascade interaction can be taken as representing the interaction of inflow distur- 
bances with the fan or the interaction of the fan wakes with the stator 6 ’ 7 . In the simplest gust-cas- 
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cade model the blades are represented as two-dimensional flat plate cascades stacked in the 
spanwise direction. The inflow distortions or the fan wakes are represented by vortical gusts con- 
verted with the mean flow. The interaction of the gust with the blades produces unsteady pressure 
distribution over the blades’ surfaces. If this unsteady surface pressure is coupled to the propaga- 
tion duct modes inside the nacelle, the radiated inlet and exhaust sound can be obtained. Classical 
work on the gust cascade problem has relied on approximate analysis and empirical formulation 
and is reviewed by Huff in Ref. 5. The two- dimension gust-cascade interaction has been studied 

by Kaji and Okazaki 8 . Three-dimensional effects are considered by Namba who pointed out that 
3-D effects become less important at high acoustic-response frequencies 9 . Effects of loading and 
nonuniform transonic flow are given by Atassi et al. 10 . Effects of camber, thickness and stagger 
angles are given in Lorence and Hall 11 . 

In the present decade considerable progress in the relatively new science of Computational 
Aero-Acoustics (CAA) has been achieved. In CAA, the unsteady governing equations are dis- 
cretized and solved for time-dependent flow variables, which includes the mean flow and the flow 
or acoustic disturbances. High-order schemes are required for discretization to reduce dissipation 
and dispersion errors. But high-order schemes support the formation of spurious modes at the 
boundaries of the computational domain. Careful attention for unsteady boundary treatment is 
needed to produce the physically correct disturbance field. 

Considerable progress has been made in applying CAA to jet noise (see review articles by 

Mankbadi 12 and Tam 13 ). Application of CAA to fan noise is still limited, and it is required to 
demonstrate first the ability of CAA to accurately predict simple models of fan noise elements. 
To achieve this purpose the gust-cascade model was considered as benchmark problem in second 

workshop on bench mark problems for CAA 14 . The amplitude of the gust was taken to be small; 
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thus, either the linearized or the full Euler equations can be used and comparison with analytical 
theories can be made. The objective was to demonstrate what CAA algorithms can offer as even- 
tual substitute for simpler source models. Several algorithms were used to compute this 

problem 15 * 17 , and the results were compared to a solution from Hall 18 . The computed results, 
however, varied considerably from Hall’s solution and from each other. 

It is suspected that the unsteady boundary treatment could be responsible for the discrepancy 
between the above simulations of the same gust-cascade problem. This might be particularly true 
for the time-domain simulations, where no exact boundary treatment exists. This issue is 
addressed in this paper. We show that the inflow and outflow treatment has considerable influence 
on the solution. Some boundary treatments are shown to produce satisfactory agreement with the- 
ory given a longer computational domain. 



Figure 1: Turbomachinery Noise Benchmark Problem 
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Problem Statement 


This benchmark problem consists of four flat-plate airfoils in an unstaggered cascade, shown 
in Fig. 1. In this problem, a periodic vortical gust is converted in a uniform mean flow, impinging 
on the flat-plate cascade and radiating noise upstream and downstream. At the inflow boundary, 
the mean flow and vortical gusts are defined as: 



C oo 


U g us, = -^cos(cu + py-G>/) (1) 

= A cos (our + Py- cot) 


where 


A = 0.01 


(X = (3 = CD = y 


( 2 ) 


The computational domain is limited to: 

-2 < x < 3 (3) 


0<y<4 



This problem is to be solved with either the linearized or the full Euler equations. The output 
requested is: (1) the real and imaginary component of the pressure jump across the airfoil, and (2) 

the intensity of the radiated sound p at the inflow and outflow boundaries. 

Governing Equations and Numerical Formulation 

The conservative 2-D Euler equations in Cartesian coordinates are given as: 
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where 


p = (Y-1)(e-^P(“ 2 + vV) ( 5 ) 

The numerical solver is a high-order accurate finite-difference time-domain scheme. The spa- 
tial derivatives are computed using a split sixth-order compact scheme 19 , and the time integration 

is performed using Hu’s 5-6 Low Dispersion and Dissipation Runge-Kutta method . The grid 
used is a uniform Cartesian grid where Ax = Ay = 0.02. Since the vortical gust has a wavelength 
of 0.8, this gives 40 points per wavelength, which produces a linear dispersion error of l.e-9. 

For all calculations, the time step used was At = 0.0125, which corresponds to a CFL number 
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of 1.2. 


Boundary Condition Formulation 

1) Tho m ps on B C 

The Thompson boundary condition 21 ' 22 uses 1-D characteristic theory to determine the 
proper boundary condition specification. In the Thompson method, we only consider the effect of 
the characteristics normal to the boundary; the tangential derivative terms are considered to be a 
source term. Thus, Eq. (4) becomes: 
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( 6 ) 


The normal flux term is then linearized and decomposed into 1-D characteristics: 


(c x ) t + (u-c)(p' x -pcu' x ) = 0 


(c 2 ), + B(c 2 p',-pV = 0 

(Cj^ + uv’, = 0 
(c 4 ) ; + (n + c)(p' x + p cu\) = 0 


( 7 ) 


where the overbar terms denote mean quantities. Thus, at the outflow boundary, 
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p' = p-p 


u' = u-u 


( 8 ) 


v' = v-v 


p' = p-p 


However, the incoming vortical gust must be specified at the inflow boundary. To accomplish 
this. 


P' = P-P 


u' = u-u-u gust 


V = v-v-v 


gust 


P' = P-P 


(9) 


These characteristics given in Eq. (7) represent, respectively, the left-running acoustic wave, 
the entropy wave, the vorticity wave, and the right-running acoustic wave. In Thompson’s 
approach, the direction of travel of the characteristic determines if it is propagating into or out of 
the domain. The characteristics that are propagating out of the domain are left unchanged; those 
propagating into the domain are set to zero. 

Once the characteristics are determined, the time derivatives of the primitive variables are 
given as: 
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(2(c 2 ), + (cj), + (c 4 ) f ) 
2c 2 

(c,), + (c 4 ), 

Pt 2 

2pc 

v r = ( c 3>, 


( 10 ) 


With the primitive variables now known, the time derivative of the conserved variables are cal- 
culated as: 


(pu), = pu, + up t 
(pv) f = pv, 

„ l._2. Pt 

E t = 2 {u ^t + — x + V uu t 


( 11 ) 


2) Giles BC 

The Giles boundary condition 23 resembles the Thompson boundary condition in that they are 
both based on 1-D characteristics. However, instead of decomposing the equations into normal 
and tangential components, Giles considered only the normal component and formulated correc- 
tions for the effect of the tangential component. Thus, while Thompson’s analysis was geared 
towards only normal and tangential 1-D waves, Giles’ analysis allows for waves with other angles 
of incidence. 
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To apply the Giles boundary condition, the interior equations are used at the boundary to com- 
pute the time derivatives. The 1-D characteristics are then defined as: 

( c i), = P't-Pw't 

(c 2 ), = c 2 p' t -p\ (12) 

(c 3 ), = P fv 'r 

(c 4 ), = P' t + pcu' t 

The outgoing characteristics are left unchanged, but the incoming characteristics are modified. 
At the inflow boundary, 

(c 2 ) y = c 2 p' y -p' y 

( C 3)y = P ?V \ ( 13 ) 

(C 4 ) y = P'y + P CU'y 

The time derivatives of the characteristics are computed using: 

(c 2 ), + v(c 2 ) y = 0 

(c 3 ) f + v(c 3 ) y + ^[(u + c)(c 4 ) y + (M-c)(c 1 ) y ] = 0 

(c 4 ), + (v)(c 4 ) y -i[(« — c)(c 3 ) y ] = 0 (14) 
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To account for the incoming vortical gust, the variables are decomposed as shown in Eq. (9). 
Once the time derivatives of the characteristics are known, Eqs. (10) and (11) are used to deter- 


mine the time derivatives of the conserved variables. 

At the outflow boundary, only one characteristic is corrected. The variables are decomposed 
as shown in Eq. (8), and the tangential derivatives are taken: 


= P'y-Pc u 'y 

(c 3 ) y = pcv’ y (15) 


The incoming characteristic is updated as: 

( C i), + “( C 3), + Kc 1 ) y = 0 (16) 


31 Hagstrom BC 

Again, like the Giles equations, the Hagstrom boundary condition 24 is formulated in terms of 
the 1-D normal characteristics: 


(c,) = p'-pcu’ 

(c 2 ) = c 2 p'-p’ (17) 


(c 3 ) = pcv' 


10 


1 I! 


(c 4 ) = p' + pcu' 


As in the Giles boundary condition, the incoming characteristics are corrected while the out- 
going characteristics are left unchanged. However, unlike the Giles boundary condition, the cor- 
rections are determined by solving additional equations at the boundary, with the user determining 
how many equations are needed for accuracy. The effect of these equations is to match to increas- 
ing order the short time expansion of the exact boundary solution. 

The inflow boundary condition takes the following form: 

( c 2 ) t = 0 

(c 3 ) / + pcjuu'j, + vv' y + ip' J = 0 (18) 

N 

(C 4 ), + p ( c 4>j, + ~ “) V ’y + \ X Vj + g j ) = 0 

i= i 

Here, f and g are the correction equations, which are given by: 

(/,), + (V + C)(fj) y + S(c t ) yy = 0 (19) 

(gj),+(y-c)(e j ) y +s(c 4 ) n , = o 

where 


it 



( 20 ) 



,_2 _ 2 . 
_ (C -M ) 

2N + 1 



For the inflow condition, the prime terms are defined in Eq. (9). 
The ouflow condition is written as: 


N 

(ci) f + v(ci) y + pc5v' y + X (fj + Sj) = 0 

7=1 


( 21 ) 


Here, the correction equations are given by: 

(/,), + (v + C)(fj) y + S(p') yy = 0 (22) 

t gj\ + (t-C)(gj) y + S(p') yy = 0 

For the outflow equations, the prime terms are defined in Eq. (8). 

4) Perfectly Matched Laver 

In the Perfectly Matched Layer (PML) boundary condition, a numerical ‘sponge layer’ is con- 
structed which surrounds the domain of interest. When a wave propagates into this layer, it is 
damped with little reflection regardless of the angle of incidence or the wavelength of the incident 
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wave. 


The PML used here is that of Abarbanel and Gottlieb 25 . In their approach, the Euler equations 
are linearized about a uniform mean flow in the x-direction. A transformation is applied to the 
equations, and a PML layer is constructed mathematically. For the cascade problem, the sponge 
layer is only in the x-direction. The PML equations are given as: 

u' t + Mu' x + p' x = -o(u' + Mp') 
v', + Mv' x + p' y = - a(v' + 2 Q + aY) + o x MP 
p' t + Mp' x + u' x + v' y = -o(Mu' + p') 

Q, + (l-M 2 )p' y = 0 (23) 

P t = (1 -M 2 )(v'-gP) 

'P, = (1 -M 2 )Q 

Here, the a denotes the damping factor, which varies through the layer as: 

a(x) = 4of V (24) 

' yX end x starl ' 

To implement these equations into the nonlinear code, the equations are transformed into: 

KW = -oW + Mp') 
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(v',) PML = - °( v ' + 2Q + a 1 ?) + g x MP 

(p\) PML = + p ') 

Q t + (1-M 2 )p’ y = 0 (25) 

P t = (1 -M 2 )(v'-oP) 

= (1 -M 2 )Q 

As before, the prime quantities are determined using Eq. (8) at the outflow, and (9) at the 
inflow. The Abarbanel and Gottlieb formulation does not include the density equation; to com- 
plete the equation set, the flow is assumed to be isentropic: 

*V.W = (P',h»L (2«) 

Once the PML corrections are calculated, they are converted into conserved variables using 
Eq. (1 1) and added to the time derivatives computed using the nonlinear Euler equations. 

In all calculations, a 20-point PML layer was employed. 

Numerical Results 

Two tests were performed for each boundary condition. In both tests, the mean flow was ini- 
tialized as: 


P = 1 
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i 1 1 1" 



u = 0.5 


v = 0 



( 27 ) 


The perturbations were initialized as: 

p g ust( x > y> °) = 0 

u gust (x, y> °) = -q(x)^cos(ax + Py - (at) 
v gust (x, y, 0) = -q(x)Acos(ax + $y - ©/) (28) 

P g usA x > y* °) = 0 


where 


q(x) = 


1 x<-2 

i + icos(j<*+2))-2<x<0 
0 x>0 


(29) 


In the first test, the cascade is removed from the domain, and the ability of the boundary con- 
dition to allow the flow to reach a periodic state is tested. In the second test, the cascade is placed 
in the flow and the result is compared with the semi-analytic result. 
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Figure 2: Pressure Time History at Inflow Boundary Point (-2,0) 
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Figure 3: Pressure Time History at Outflow Boundary Point (3,0) 
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0 Test Problem 1: No Cascade 


In the first test problem, the domain was initialized with the perturbation given in Eq. (28), 
and the cascade was removed. This tests the implementation of the upstream disturbance as well 
as the ability of the outflow boundary condition to allow the vortical gust to convert out without 
reflection. The grid for all cases had a uniform spacing of Ax = Ay = 0.02. 

In these problems, the time history of pressure is shown at two points, (-2,0) and (3,0). Figures 
2 and 3 show the pressure time history at these two points for the different boundary conditions. 
From top to bottom, the graphs show the Thompson BC, the Giles BC, the Hagstrom BC using 3 
equations, the Hagstrom BC using 10 equations, and the Perfectly Matched Layer (PML) BC. 
Notice that the scales are all similar except for the Thompson boundary condition, which does not 
give an adequate solution for this problem. Also, the 10-equation Hagstrom condition does not 
give a better solution than the 3-equation Hagstrom boundary condition. This is not surprising, as 
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the Hagstrom BC is designed to perform well for transients and not for long-term time marching, 
and, as the number of auxiliary equations is increased, the long-term solution may be less accu- 
rate. 

Figures 4 and 5 show the v velocity and pressure distributions along the y=0 line. The inflow 
vorticity specification is working well, as shown in Fig. 4. Figure 5 illustrates the pressure error, 
which is mainly due to the boundary conditions reflecting the outgoing vorticity and the acoustic 
waves generated by the initial vorticity specification. From Figure 5, the Thompson condition is 
giving a much worse solution than the other boundary conditions; the other boundary conditions 
have roughly equivalent levels of error. 


21 Test Problem 2: Specified Domain Calculation 

The next calculations were performed on the specified domain, which covered -2 < x < 3. All 



Figure 5: Pressure Distribution along y = 0 Line 
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of the boundary conditions were tested for this configuration. 

Figures 6 and 7 show the pressure time history at the inflow and outflow boundary points. 
Some of the data is fragmentary; this is due to the lead author’s inability to type while transferring 
data. 

The Thompson boundary condition proved to be unstable due to reflections from the inflow 
boundary interacting with the cascade and initiating a feedback loop. Figure 6 illustrates the rapid 
buildup of pressure in the domain. 

The Hagstrom boundary conditions were very slowly unstable for this problem, as Figs. 6 and 
7 show. As expected, adding more equations at the boundary improved the transient solution; 
however, the long-term periodic solution was not improved. Also, as more equations were added, 
the solutions took longer to converge. Again, it must be noted that the Hagstrom boundary condi- 
tion was developed for transient problems and not for long-term periodic problems. 



Figure 6: Pressure Time History at Inflow Boundary Point (-2,0) 
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Both the Giles and the Perfectly Matched Layer boundary conditions allowed a stable long- 
term solution to be calculated. The PML condition required extra points in the buffer zone; how- 
ever, it did improve the solution at the boundary somewhat. 

Figures 8 and 9 show the RMS pressure results calculated using these boundary conditions. 
All solutions were in the same range, with the PML boundary condition performing slightly better 
than the rest. 

Figure 10 shows the pressure difference across the y=0 airfoil for the Giles boundary condi- 
tion; all of the stable boundary conditions allowed identical airfoil results to be obtained. The 
effect of the numerical singularities at the leading and trailing edge are illustrated. Wall conditions 
were applied for the upper and lower leading and trailing edge points, with the results being aver- 
aged' and applied to both points. 

It must be noted that these singularities are purely due to the infinitely thin airfoils prescribed 
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Figure 7: Pressure Time History at Outflow Boundary Point (3,0) 
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Figure 9: RMS pressure at downstream boundary 
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in the problem statement; it is expected that these singularities are the main factor causing the 
error on the airfoil. Future work will address the effects of real airfoil geometries. 

3^ Test Problem 3: Long Domain Calculation 

To quantify the relative effect of the numerical singularites and the close boundaries on the 
computed solution, the same calculation was performed on a long domain (-8 < x < 9) using the 
same grid spacing and time step. For these calculations, which required much more CPU time to 
converge, only the Giles and PML boundary conditions were used. Figure 1 1 shows that the solu- 
tion on the sirfoil does not change as compared to the short domain solution. 

Figures 12 and 13 show the computed RMS pressure disturbances at the upstream and down- 
stream locations. Using this long domain, both boundary conditions give a solution that is close to 
the exact solution in phase, with a greater amplitude. It is theorized that the errors in the pressure 



Figure 10: Comparison of Airfoil Pressure Distribution with Exact Solution 
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on the airfoil due to the numerical singularities are causing the amplitude error. However, the 
PML outflow condition is performing better than the Giles outflow condition, as shown in Fig. 13. 

Figure 14 shows the effect of the extended domain on the pressure contours using the Giles 
boundary condition. On the airfoil itself, the contours are nearly identical. In the upstream direc- 
tion, the contours become different as the inflow boundary is approached, showing that there is 
some reflection from the inflow boundary condition. In the downstream direction, however, the 
contours diverge almost immediately after the trailing edge, suggesting that the outflow boundary 
condition is more reflective than the inflow. 

Conclusions 

Several boundary conditions were tested on a CAA Workshop Benchmark Cascade problem. 
The problem is very difficult due to numerical singularities in the problem specification as well as 



Figure 11: Comparison of Airfoil Pressure Distribution for Long Domain Solution 
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the proximity of the boundaries. Since the problem in periodic in one direction, the ability of the 
boundary condition to allow acoustic and vortical waves to exit without reflection determines the 
stability and accuracy of the solution. It was found that the error due to the numerical singularities 
at the leading and trailing edges of the flat plates kept the solver from obtaining the exact ampli- 
tude of the radiated sound; however, the reflected waves at the boundaries due to the boundary 
conditions kept the solver from obtaining the correct radiation pattern even though the airfoil 
pressure distribution was nearly identical. 

Of the four boundary conditions tested, the Thompson condition was found to have the high- 
est reflection at the boundaries, resulting in code instability. The Giles condition did reasonably 
well while requiring very little extra computational work. The Hagstrom condition did not per- 
form well for this problem, but it must be noted that it is designed for short-term transient calcula- 
tions and not long-term periodic problems such as this one. The Abarbanel and Gottlieb PML 



Figure 14: Comparison of Specified-Domain and Long-Domain Instantaneous Pressure Contours 

with Giles Boundary Conditions 
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boundary condition performed the best, but required a 20-point sponge layer at each end. How- 
ever, the performance of the PML condition could probably be improved with better specification 
of the damping coefficient distribution in the sponge layer. 

Currently, boundary condition specification is one of the main sources of error in time-domain 
CAA calculations. Work in this area will yield large dividends. 
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